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Abstract: We analyse a specific, duality-based generalization of the hermitean matrix 
model. The existence of two collective fields enables us to describe specific excitations of 
the hermitean matrix model. By using these two fields, we construct topologically non- 
trivial solutions (BPS solitons) of the model. We find the low-energy spectrum of quantum 
fluctuations around the uniform solution. Furthermore, we construct the wave functional 
of the ground state and obtain the corresponding Green function. 
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1. Introduction 

The recently renewed interest in matrix models came from the identification ^ of ma- 
trix quantum mechanics with the worldvolume theory of unstable DO-branes in the two- 
dimensional string theory. One uses matrix models to study non-perturbative aspects of 
string theory, or more specifically, the open-closed string duality. 

In the usual approach, the closed string states are related to the matrix model states by 
bosonization of the fermion field followed by a non-local field redefinition Q . Alternatively, 
one can use the collective- field formulation of the matrix model to obtain a continuous, field- 
theoretic description of the model. This approach was used to construct the collective-field 
formulation of unitary Q and symplectic j|] matrix models. The important result of the 
collective-field formulation was the analysis of higher-order terms in the 1/iV expansion. 
It was shown that divergences cancelled and the finite correction, of order 1/N 2 , was 
calculated || . 

The two-dimensional string theory is obtained from a matrix model in the double 
scaling limit, with the critical potential given by the inverted harmonic oscillator. It 
was shown in Ref.||] that one could remove the oscillator potential from the action by 
a suitable coordinate transformation. Therefore, it is of interest to study general features 
of the background independent theory without the external potential. 

Furthermore, it was realized that in order to describe black-hole solution, the general- 
ization of the standard hermitean matrix model is needed Q . One such generalization was 
proposed in the paper []|] . 

In this paper, we analyze this duality-based generalization |J] and interpret the model 
in question as hermitean matrix model beyond the usual singlet approximation. We con- 
struct non-trivial solitonic solutions of the model, and find the spectrum of quantum fluc- 
tuations around the uniform density solution of the master Hamiltonian constructed in 
Refs.|,|]. 
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The master Hamiltonian can be interpreted as a two-family Calogero model, connected 
by duality, with mutually inverse coupling constants J?]], or, for a special choice of coupling 
strength (A = 1/2), as a master Hamiltonian for symmetric and quaternionic models ||. In 
Ref.|| it was conjectured that the master Hamiltonian with A = 1/2 actually corresponded 
to the hermitean matrix model, and our analysis of the model confirms this conjecture. 
From the string point of view, we describe oriented string degrees of freedom (hermitean 
matrix model) by unoriented string degrees of freedom (symmetric and quaternionic matrix 
model) 1 . Although we focus on the matrix model with A = 1/2, constructed in Ref.||, 
throughout the paper we write A for generality. Also, in this way, the A <-> 1/A duality is 
obvious. The aim of this investigation is to shed more light on the nature of duality in the 
field-theoretic description of matrix models. 



2. Semiclassical solutions of the master Hamiltonian 



Let us start with the master Hamiltonian obtained in Ref.fl 



if " = i / dxp(x)(d x * f (x)) 2 + i / ,/..,.[..■) 



p{x) 



2\hy^ + 2-1 dy 
J x-y J 
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The collective fields p(x) and m(x) are normalized as 



(2.1) 



dxp(x) = N, / dxm(x) = M, 



(2.2) 



where N and M are large numbers of independent matrix eigenvalues (or particles in the 
system for the Calogero model), and vr p / m (x) represent the canonical momenta: 



[d x TC p (x), p(y)] = -id x 5(x -y), [d x ir m (x),m(y)] = -id x 5(x -y). 



(2.3) 



One can add a term —p p J dxp(x) and an analogous one for m(x) to the Hamiltonian to 
ensure Q2.2|) . 



The leading part 2 of the collective-field Hamiltonian in the 1 /N and 1 /M expansions 
is given by the positive-definite effective potential, the second and the fourth term in the 



Hamiltonian (2.1). Owing to the positive definiteness of the effective potential, the Bogo- 
mol'nyi limit appears. The Bogomol'nyi bound is saturated by the positive normalizable 



1 One can think of our master Hamiltonian as an explicit realization of an idea proposed in pX] , that 
one can continuously deform oriented string theories in two dimensions into unoriented ones by turning on 
non-local interactions on the worldsheet. 

2 Note that the divergent terms (the last line in the Hamiltonian (2.1), as well as the momenta, are 



suppressed in the leading order, but are contributing to the spectrum of quantum fluctuations. 
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solutions po(x) and mo(x) of the coupled equations 

(A - + 2\JdvM. + = 0, (2.4) 

p(x) J x - y J x-y 

n_\ dMx) + u + 2 r dy Ay±_ = (2 5) 

\A J m[x) XJ x-y J x-y 
An obvious solution is a constant density for both fields: 

po = N/L, and m = M/L. (2.6) 

This solution exists on the compact support only, as can be seen from the normalization 
conditions (|2.2j). 

Next, let us look for non-trivial solutions of the coupled equations ( p.4[ ) and (|2.5|). 
Substracting these equations and then integrating them, we find the condition 

p{x)m{x) = c, (2-7) 

where c is an arbitrary constant parameter. Expressing m(x) and introducing it in eq. 
dU), we are left with 



(A - l)a^(x) + 2\p(x)-l d ^- + 2cp(x)-f—-pt = 0. 

J x-y J p{y)(x - y) 



(2.8) 



Let us take the following ansatz 



x 2 + a 2 



x z + o z 



where a and b are arbitrary positive constants. Using the Hilbert transform 



x — y y — z x — z 

we find the conditions 



dy 1 siqnilm z) 

iTT^-± '-, (2.10) 



b (^-^)= V' c= ^- (2 - n) 

The soliton-antisoliton solution is given by 

p(l)=w + i_I_JL_, M 

, . c 1 — A o , 

m(x) = — + 5 5 . 2.13 

po tt x A + or 

Depending on the sign of (A — 1), one density describes a hole in the condensate, whereas 
the other describes a particle above the condensate. We see that A <-> 1/A interchanges 
holes and particles. 

Now let us study the case c = ka in the limit a —* 0. In this limit, m(x) goes to the 
delta function: 

lim m(x) = 7T(5(x) (2-14) 

a^o p 
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and p(x) reduces to 



limp(x) = p 

a—>0 X z + Z 



(2.15) 



Although these two expressions are obviously solutions of the Bogomol'nyi equations, it is 
instructive to insert them into eqs. (|2,4|), ( |2.5| ) to confirm the correctness of the limiting 
procedure. Inserting ( 2.14| ) and ( 2.15) ) into eq. ( |2.4j ), we obtain the conditions 



k = ^ , irp b 



1 - A 
A ' 



Next, we check the validity of eq. 



dy h 2Xm(x)4- dy 

x-y J 



p(y) 

x-y 



Inserting ([D3]) and pd5| ) into ( pH) ) we find 



ka 2x(x 2 + a 2 ) - 2x(x 2 + b 2 ) £;a x 2 + 6 2 7rp (a 2 - b 2 ) 
(1 — A) lim ; n — — h 2A lim - - hm 



+2 lim 



a->0 p 
fca x 2 + b 2 



lim 



(x 2 + a 

7Tfc(6 2 



2^2 



(2.16) 



(2.17) 



+ 



a-»0 po x 2 + a 2 a^O 

Using the identity 



po £ 2 + a 2 



lim ■ 



2to P 

2 = 27T — 0(X) 



( x 2 + 6 2) 

and ( |2.14| ), we obtain only one condition 

irkb 2 
Po 



x 



-7rd x 5(x) 



1- A 



x 2 + b 2 
(2.18) 

(2.19) 
(2.20) 



which is obviously consistent with ( [2.16 ). This is due to m(x) ~ <5(x) and ( |2.17| ). 



Next, we construct the two-soliton solution of eqs. (|2.4[) , (2.5). With the following 
ansatz for p(x) and m(x) 

(x — a)(x — a*)(x + a)(x + a*) 



p(x) 



m(x) 



(x-b)(x-b*){x + b)(x + b*) ' 



m - 



(x - 6)(x - 6*)(x + 6)(x + b*) 



(2.21) 



p (x — a)(x — a*)(x + a)(x + a*) ' 

we easily find the condition b = a, so that the solutions reduce to trivial ones with constant 
po and mo. 

Now, let us take c = klm(a) in the limit Im(a) — ► 0, and write the ansatz ( 2.2 1] ) in 
the more natural form 



p(x) 



lim po 

e— >0 



[(x - x ) 2 + e 2 ][(x + x ) 2 + e 2 ] (x - x ) 2 (x + x Q f 



(x 2 — 6 2 )(x 2 



Po 



(x 2 - 6 2 )(x 2 - b* 2 ) 



B B* 
Po 1 + r + 



B 



B* 



x + b x + b* 



(2.22) 



m x 



lim 



he 



(x 2 -6 2 )(x 2 



e^o po [(x - x ) 2 + e 2 ] [(x + x ) 2 + e 2 ] 



k [S(x - x ) + 5(x + x )] , (2.23) 
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where 



B 

k 



x, 



2\2 



26(62 _ 6 *2) . 

kirjxl - b 2 ){xl - b* 2 ) 
Ap xl 



(2.24) 



Inserting the ansatz ( 2.22j ) into eq. (|2.4j ) we find the following three conditions: 

B = -B* , 2\p iriB = A - 1 , k = 1 - A. 
After we introduce b = \b\e %< ^ and r = Xq/|6| 2 > 0, the condition B = —B* gives 



r = 2 cos <j) — 1 xq = \b\ y (2 cos — 1), 



(2.25) 



(2.26) 



Next, from pqttiB = A — 1 we find 



p vr|6| 



1 — A cot | 
A 2 cos d> 



(2.27) 



So, from ( |2.25 ) and ( 2.26| ) we have obtained \b\ and xq as functions of <fi and the last 
condition gives k 



2klT 

Po 



cos 



(1 — COS < 



2 cos i 



1 



1-A, 



(2.28) 



as a function of 4>. 

In the BPS limit, one can construct n-soliton solutions, n > 2, following the strategy 
demostrated above for the two-soliton solution. The general ansatz for n-soliton solution 
is of the form 



i=l 



P (x)= P Ai + Y,{ — + B ' 



a=l 
rn 



rf >y rp ryf 

GL % * j O 

2 



(2.29) 



Note that this ansatz satisfies the condition p{x)m{x) = 0. To obtain relations between 
poles and zeros explicitly, one has to solve the general algebraic problem of finding zeros of 
the polynomial of degree n > 3. The soliton solutions have the following properties: they 
are located at the point at which the p(x) field is vanishing, but the m{x) field is becoming 
undetermined. This property is also characteristic of BPS solutions describing monopols 
and Julia-Zee dyons [11]. 
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3. Quantum excitations 

3.1 The spectrum of quantum fluctuations 

At this point we analyse the dynamics of the collective-field excitations around the ground- 
state solution (|2.6| ) of our master Hamiltonian fl2.1|) . First, we rewrite the master Hamil- 
tonian in the following form: 



H 1 ' 



M = l [ dxp(x)Al(x)A (x) + £ ( dxm(x)Bl(x)B (x), 



(3.1) 



where 



A) 0*0 = d x -K p (x) + i 
B Q (x) = d x n m (x) + i 



(A - 1) d xP (x) 

2 p(x) 
(1 — A) d x m(x) 1 f m(y) 



J x-y J x -y 

p(y) 



+ 



dy- 



+ +dy 



x-y 



2A m(x) Xj x — y 
Next, we perform the 1/N (1/M) expansion of the collective field p(x) {m(x)) 

p(x) = po + d x r](x), m(x) = m + d x rj(x), 



(3.2) 



(3.3) 



where d x rj{x) and d x fj{x) are small density quantum fluctations 3 . We insert (|3.3[) into the 



Hamiltonian (3.1) and expand up to second-order terms in d x r/(x) and d x rj(x). We obtain 
a Hamiltonian quadratic in fluctuations: 



H m = PJL j dxA \ A + dxB^B. 
The operators A(x) and B{x) 



(3.4) 



A(x) 
B(x) 



-7T v (x) + i 
-TTrj(x) + i 



(A-l) ^r? 
2 /J 

2A mo 

satisfy the following commutation relations: 

A-l 



+ 



+ 



dy 
x-y 

dy 



iXd y rj(y) + d v fj(y)) 

_ {dyv(y) + \d y fj{y)) 
x y a 



A{x),A\y) 
B(x),BHy) 
A{x),B\y) 



d x d y 5(x -y) + 2Xd x 



P 



20 



Po x-y 

d x d y S(x -y) + -d x - 

Xitiq 

P 



A x — y 



x-y 



(3.5) 
(3.6) 

(3.7) 
(3.8) 
(3.9) 



Now, to find the spectrum of low-lying excitations, we have to diagonalize the Hamil- 
tonian ( |3.4|) . We expand the operators A(x) and B(x) in terms of new, complete sets of 



3 The small density quantum fluctations are defined by the explicit derivative because of the normalization 
condition J dxd x rj(x) = 0. 
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operators 



n n 



where 



0"ni O-ra 



LO n b r , 



b 



On; 



bt 



^0. 



We demand that the Hamiltonian ( |3.4D should take the following form: 

= Y^aia n + Y,bib n . 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



Generally, n is the quantum number, and it is assumed that the sum is replaced by an 
integral for a continuous spectrum. 

Let us concentrate on the first part of the Hamiltonian ( |3.13D . We insert the expansion 
( 3.10 ) in the commutators ( |3.7| ) and apply ( [3.12 ) to obtain the completeness relation 



^w n ^„(x)0*(y) 



A — 1 P 
-d x d y S(x -y) + 2Xd x 



Po x-y 
Inserting (3.10) in ( |3.4| ) and demanding ( p.13 ), we obtain an orthogonality relation 

' dxcj)* n (x)(f) m (x) = 5 nm . 



(3.14) 



(3.15) 



Next, we multiply the relation ( 3.15| ) by ^(y), sum over m, apply the completeness relation 

(y) 



( 3.14 ) and finally obtain the equation for the functions 

A — 1 



[x)\ 



Vn4> n (x) 



-d x <p n (x) + Xp d x fdy 



x-y 



The plane waves are the solution of Eq. ( |3.16| ) with 

-n 2 + \poir\n\ 



1-A 2 



An analogous procedure for the second part of the Hamiltonian ( |3.13| ) gives us Q n 

A - 1 



-n +moir\n\. 



(3.16) 



(3.17) 



(3.18) 



A few comments are in order here. The relation (3.17) was obtained in Ref.|12] as 
a dispersion relation for quantum fluctuations around the ground-state solution in the 
Calogero model. Under the duality transformation A — > 1/A and po *~~* m o ; the relations 



( 3.17 ) and ( 3.1 8| ) are interchanged, reflecting the known duality symmetry A — > 1/A of the 
Calogero model. However, note that these dispersion relations do not describe physical 
excitations since [a, b>] ^ 0, i.e., the Hamiltonian ( |3.13 ) is not diagonal! 
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Furthermore, for A = 1/2, the relations ( |3.17| ) and ( 3.18 ) were obtained in Ref . |l3| . 
In this paper, these dispersions represented quantum fluctuations for two independent, 
effective Hamiltonians operating on two particular subsets of eigenstates of the hermitean 
matrix model 4 . However, our duality-based generalization of the hermitean matrix model 
also includes the interaction between these two sets of mutually dual degrees of freedom. 

Let us continue with our diagonalization procedure. From the completeness and or- 
thogonality relations we obtain 



exp(mx) exp(mx) 

(f>n{x) = ; , ifn[x) - 



Now, we are in a position to calculate the commutator 
using (PI) , ( pit ), and (p9[) , we obtain 



(3.19) 

From Q and (p^ ), 



n,b m = 5 nm \n\iT\f p m \ = 5 nm f n . 



(3.20) 



In the limit of small n, an interesting possibility arises - a n oc b n . Appropriate rescaling 
of the relations (|3.12j ) gives us the following Hamiltonian: 



H (2) = ^( Wn + n n )clc n = |n|7r(Ap + m )cic n , 



(3.21) 



where the bosonic operators c n satisfy 



c f c f 1 



C-rii C„ 



(3.22) 



We see that in the lowest order, the special case A = 1/2 reproduces the part of singlet 



sector of hermitean matrix model [12], i.e., we can represent the hermitean matrix model 
(the Calogero model with A = 1) as a dual system of symmetric (A = 1/2) and quaternionic 
(A = 2) matrices. 

Of course, this is to be expected. Although we claim that the collective- field variables 
used in the construction || of the master Hamiltonian ( |3.1| ) with A = 1/2 do contain 
some non-singlet states of the hermitean matrix model, these non-singlet states do not 
contribute in the low-energy approximation, and, therefore, the obtained result is in fact 
a consistency check. In the string theory language, this means that the unoriented string 
amplitudes at tree level are the same as in the oriented model. 

What can we say about the excitations for "not-so-low" energy? Let us exactly diag- 
onalize the Hamiltonian ( |3.13| ). Having the quadratic Hamiltonian allows us to apply the 
Bogoliubov transformation. We introduce two sets of mutually commuting operators c n 
and d„ 



(3.23) 



4 These effective Hamiltonians were actually a collective-field formulation of matrix models for the sym- 
metric and quaternionic matrices, compare Ref.bl. 
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such that we can write the Hamiltonian ( 3.13] ) as a sum of two independent quadratic 
parts: 

H^=J2c% + J2^dn. (3-24) 



Assume that 



- p n d n and b n = (3 n c n + a n d n , 



(3.25) 



where a n ,(5 n are some coefficients that depend on the quantum number n. Now, inserting 
(3.25) into Q3.13 ) and demanding ( |3.24j) , and inserting ( |3.25 ) into ( p. 12 ), gives us a system 
of four equations 

a n /3 n (w+ -oj~) = f n , c? n + Pi 



a 2 n uj+ + (3 2 n uj n = oj n , f3 2 u+ + . 



= 1, 



(3.26) 



which we solve. The coefficients in the Bogoliubov transformation ( 3.25 ) are 



i + 



ft = o 1 



vk-o„)2+4/2; ' 



LO,, 



VK-^„) 2 + 4/^ 
given 

+ Q n ± v/K - ^n) 2 + 4/2 



(3.27) 



The spectrum of low-lying excitations is given by 
-u 1 



7r(Ap + m )\n\ 



1 / 2(1 A)(A w - r o) 
V vr(Ap + m) 2 



(1-A)s 



7T (Apo + m ) 



r n 2 (3.28) 



Note that each branch of the dispersion relation ( 3.28 ) is invariant under the duality trans- 
formation A — > 1/A and po <-> mo, and physical excitations are not connected by duality. 
We have two sets of independent oscillators. Also note that obtained corrections are of the 
same order as non-singlet corrections [15|, namely 1/N 2 with respect to the leading singlet 
contribution. 

However, for A = 1/2 the ground system is dual, i.e., Xpo = mo and uj~ is negative! 
The "wrong" sign of iv^ signalizes the instability of the vacuum. 



3.2 The vacuum functional and Green functions 

The vacuum functional is a simultaneous solution of the equations 

Aq(x)\0) = and B (x)\0) = 0, 

and is given by 



(3.29) 



*o(p,™) = exp 



+ 



1) 



dxp{x) In p{x) + — 



dxdyp(x) In \x — y\m(y) + 



(1-A) 
2A 



dxm(x) lnm(y) + 



2A 



dxdyp(x) In \x — y\p{y) + 



dxdym(x) In \x — y\m(y) > . (3.30) 
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We expand the vacuum found in ( |3.30 ) up to quadratic terms 



*( 2 ) = exp (J j dzdy rj T (z)G- 1 (z - y)r,(y)^ , 



(3.31) 



where we have introduced the matrix notation 



(l—\)nS(z—y)—2\irpo In \z— y\ — 2iv^J\pQmQ In \z— y\ 

— 2TTy/\pomQ In \z— y\ (X—l)irS(z—y) — 2iTm() In \z— y\ 



Jk{z-y) 



~& / dk — 



^>k fk 
fk ^k 



This vacuum is a solution of the equations 

A(x)^ 2) = 0, B(x)¥ 2) =0« c fc ^ (2) = 0, 4^ {2) = 
Writing in terms of the field 



(3.32) 



(3.33) 



<p(x) = J dxS(x — x')ri(x'), 
S(x — x') = 



ik(x~x') 



dk- 



27r^/\k\(k 2 + A; 2 ) 1 /* 
where /cq = 4A 2 7r 2 / o 2 l /(l - A) 2 , we obtain 



Ok Pk 

-Pk ot k 



(3.34) 



exp 



dzd yv > T (z) (-— [ dfce ifc (-y)V^+M 

V 4?r J \ k \ 



o u7 



<P(V)> (3-35) 



and the non-normalizability of the vacuum vl^ 2 ) is manifest. This is a consequence of 



the commutation relations ( p. 23 ) from which we conclude that d n is actually a creation 
operator which we call dn and d n is an annihilation operator now called d n . 
The Hamiltonian ( |3,24| ) is, after appropriate rescaling, given by 

ij( 2 ) = ^u+f^cn + ^ l w n \ d n4i = ^W+C+C n + ^ \ d l d n + ^2\^n\- (3-36) 

ran n ra ra 

The last term must be included in the vacuum energy. It simply defines the Fermi level of 
the system. 

Solving equations for the normalizable vacuum functional 



c k $M = ,4$ (2) = 0, 



we find 



$( 2 )=exp|| J dzd yi p T {z)G-l{z-y) V {y)Y 



(3.37) 



(3.38) 
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The inverse propagator G^(z — y) is given by 



1 

4 

A-l 



G'hz-y) = -- I dke* k ^ 



K7| 



(dzd y + kl) 



(l*)5{z-y) + (l\) k -±K {k,{z-y)) 



, (3.39) 



where Kq{x) is the Bessel function. 

The two components of the field <p(y) can be combined into a complex scalar field, 
which could be interpreted as a complex tachyon field of 0^4 two-dimensional string the- 
ory ||]. 



4. Conclusion 

We have analysed the spectrum of quantum fluctuations around a particular ground-state 
solution of the master Hamiltonian Q3.1| ) , representing a duality-based generalization of the 
hermitean matrix model. This has given us a glimpse of the dynamics of dual degrees of 
freedom in field theory. We have shown that for A = 1/2, the master Hamiltonian describes 
a new representation of the hermitean matrix model, in which some specific states have 
been analysed. These states contribute to the next-to-leading order, and the contribution 
is comparable with the contribution from the non-singlet sector. Also, we are able to 
construct topologically non-trivial, solitonic solutions of the model in question, which do 



not exist in the singlet approximation of the hermitean matrix model. In Ref.[15| the 
non-singlet states were attributed to the vortices on the worldsheet. Note that these non- 
local interactions on the worldsheet can be interpreted as double-trace deformations of the 
standard hermitean matrix model originating from the (A — l)-proportional terms in the 



master Hamiltonian [10|. 

The excitations in our model behave as quasi-particles and quasi-holes, and can be 
described by a complex scalar field. Therefore, in our opinion, the relation between the 
matrix model constructed in Ref.Q and OA string theory in two dimensions deserves a 
more detailed analysis. Furthermore, we have constructed the soliton-antisoliton solution 
for the master Hamiltonian, and quantum fluctuations around this solution are currently 



analysed [19]. Also, it would be interesting to investigate the possible integrability of the 
model, as the existence of multisoliton solutions suggests. 
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